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Non-baryonic, or "dark" , matter is believed to be a major component of the total mass 
budget of the universe. We review the candidates for particle dark matter and discuss 
the prospects for direct detection (via interaction of dark matter particles with labora- 
tory detectors) and indirect detection (via observations of the products of dark matter 
self-annihilations), focusing in particular on the Galactic center, which is among the 
most promising targets for indirect detection studies. The gravitational potential at the 
Galactic center is dominated by stars and by the supermassive black hole, and the dark 
matter distribution is expected to evolve on sub-parsec scales due to interaction with 
these components. We discuss the dominant interaction mechanisms and show how they 
can be used to rule out certain extreme models for the dark matter distribution, thus 
increasing the information that can be gleaned from indirect detection searches. 

1. Dark Matter Candidates and Constraints 

If one compares the total amount of matter inferred from cosmic microwave back- 
ground (CMB) experiments^ 

fl M h 2 = 0.135™. (1) 

(here £1 denotes the mean density as a fraction of the critical, or closure, density, 
and h is the Hubble constant in units of 100 km s _1 Mpc -1 ) with the amount of 
baryons allowed by the CMB 

fl b h 2 = 0.0224 ± 0.0009 (2) 

(also consistent with the Big Bang nucleosynthesis constraint 0.018 < Vt^h 2 < 0.023, 
see e.g. Ref. |2J), one is left with a component of the Universe that is composed of 
matter distinct from ordinary baryonic matter. There is no shortage of dark matter 
candidates, most of them arising in theories beyond the standard model of particle 
physics. In the framework of the standard model, neutrinos have been proposed 
as dark matter candidates, but the analysis of CMB anisotropies, combined with 
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large-scale structure data, suggest Sl v h 2 < 0.0067 (95% confidence limit), which 
implies that neutrinos are a sub-dominant component of non-baryonic dark matter. 

Rather than compile a complete list of all possible dark matter candidates (often 
referred to as WIMPs: "Weakly Interacting Massive Particles" ) we discuss here the 
candidates that have received the widest attention in the recent literature. For an 
extensive review of particle dark matter candidates see e.g. Refs. EUand^J 

The neutralino in models of R-parity-conserving supersymmetry is by far the 
most widely studied dark matter candidate. In the "minimal supersymmetric stan- 
dard model" (MSSM), the superpartners of the B and W3 gauge bosons (or the 
photon and Z, equivalently) and the neutral Higgs bosons, H® and H® , are called 
binos (B), winos (W3), and higgsinos (H® and H^), respectively. These states mix 
into four Majorana fermionic mass eigenstates, called neutralinos. The lightest of 
the four neutralinos, commonly referred to as "the" neutralino, is in most supersym- 
metric scenarios the lightest supersymmetric particle. The neutralino is a perfect 
dark matter candidate, with mass and interaction cross sections that, in fine-tuned 
regions of the supersymmetric parameter space, can correctly reproduce the dark 
matter relic density while still being consistent with constraints on the WIMP mass 
and searches in accelerators as well as direct and indirect detection experiments. 

Another interesting candidate arises in theories with universal extra dimensions 
(UED)jSl i.e. extra-dimension scenarios in which all fields are allowed to propagate 
in the bulk. Upon compactification of the extra dimensions, all of the fields propa- 
gating in the bulk have their momentum quantized in units of p 2 ~ 1 / R 2 , where R is 
the compactification radius, appearing as a tower of states with masses m n = n/R, 
where n labels the mode number. Each of these new states contains the same quan- 
tum numbers (charge, color, etc.) The lightest Kaluza-Klein particle (LKP) in the 
framework of UED, which is likely to be associated with the first KK excitation 
of the photon, or more precisely the first KK excitation of the hypercharge gauge 
bosonpl provides a viable dark matter candidate. We will refer to this state as 
The calculation of the B^ relic density was performed by Servant and TaitPwho 
found that if the LKP is to account for the observed quantity of dark matter, its 
mass (which is inversely proportional to the compactification radius R) should lie 
in the range 400 to 1200 GeV, well above any current experimental constraint. 

Dark matter candidates are commonly believed to have masses i n th e range 1 
GeV - 100 TeV. The lower value is the so-called Lee- Weinb erg limit while the 
upper value comes from the so-called unitarity 

boundHnini 

It is sometimes stated 

in the literature that the upper bound is of order 340 TeV. However, this comes 
from the now-obsolete requirement that VtMh 2 < 1, adopted in 1989 by Griest and 
KamionkowskiPl If one instead adopts the modern estimate in Eq. (Q, the upper 
limit on the mass of the dark matter particle is m < 120TeV. 

These limits are, however, mo del- dependent. L ight sc alar particles with masses 
below 1 GeV are viable dark matter candidates 1121121 (referred to as light dark 
matter or alternatively MeV dark matter), since the Lee- Weinberg limit strictly 
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applies only to fermionic particles with standard model interactions. If MeV dark 
matter is to explain t he 511 keV emission from the Galactic bulge observed by 
INTEGRAL, l^^l^vUO l then a comparison of the inverse Bremsstrahlung emission 
(associated with dark matter annihilations into electron-positron pairs) with the 
diffuse Galactic background constrains the mass of the light dark matter particle to 
be smaller than ~ 20 MeV.^ht has been argued that hints of this scenario may also 

l-i 61 

have been discovered in particle physics experiments and could have interesting 
implications for neutrino physicsP^ 

The unitarity bound can also be violated if, e.g., dark matter particles were 
not in thermal equilibrium in the early Universe, but were in stead produced via 
alternative mechanisms, such as gravitational production! ^ A 

2. Halo Models 

The probability of direct detection is proportional to the dark matter density in the 
solar neighborhood, and most discussions of indirect detection have also focussed on 
the Milky Way; in both cases, a model of the Milky Way's dark matter halo is crucial 
for predicting and interpreting event rates. Halo density profiles are usually derived 
from A-body simulations of gravitational clustering, which predict a characteristic 
form for p{r), and on dynamical constraints like the Sun's orbital velocity, which 
provide a normalization. (For revi ews, se e Refs. |221 ~ E31) The most recent and 
highest-quality A-body simulations^*^! suggest a universal dark-matter density 
profile of the form 

/o(0 « /° exp ~(r/r ) 1/n (3) 

with n ps 5;-^ the density normalization at the Solar circle implied by the Galactic 
rotation curve is p ~ 0.3 GeV cm' 3 » 8 x 10~ 3 M©pc~ 3 , with an uncertainty of 
a factor ~ 2. Henceforth we refer to Eq. J3J) with p(Rq) = 0.3GeVcm -3 as the 
standard halo model (SHM). (Eq. replaces the more approximate NFW modeP^ 
in which the density is a single power-law, p ~ r _1 , inward of the Solar radius Rq.) 
Unfortunately the iV-body simulations from which these halo models are derived 
have resolutions that are measured in hundreds of parsecs at best, whereas the signal 
from dark matter annihilations depends critically on the mass profile in the inner 
few parsecs of the Galaxy. Furthermore the A-body simulations typically ignore the 
influence of the baryons (stars, gas, etc.) even though these are known to dominate 
the gravitational potential of the inner Galaxy. 

One simple, and probably highly idealized, way to account for the effect of the 
baryons on the dark matter is via adiabatic contraction models, which posit that 
the baryons contracted quasi-statically and symmetrically within the pre-existing 
dark matter halo, pulling in the dark matter and increasing its density - When 
applied to a dark matter halo with the density law © , the result is a halo pr ofile with 
p ~ r ,c , 7 C w 1.5 inward of the Solar circle and an increased density at i? 
Alternatively, strong departures from spherical symmetry - for instance, during the 
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mergers that created the Galactic bulge - might have resulted in lower dark matter 
densities in the inner tens of parsecsP^^2Hl However we argue below that the 
steeply-rising stellar density near the Galactic center makes such models unlikely, 
at least in the case of the Milky Way. 

The annihilation signal from a region of volume V that includes the Galactic 
center is proportional to (p 2 )V. For any p oc r -7 with 7 > 3/2, the small-radius 
dependence implies a divergent flux; hence the strength of an annihilation signal 
is crucially dependent on the dark matter distribution on very small scales, where 
neither the simulations nor the dynamical data provide useful constraints. In the 
following section we discuss physical arguments that can be used to constrain the 
dark matter distribution on these small scales, and in the final section we discuss 
the implications for indirect detection. 

3. Dark-Matter Dynamics on Sub-Parsec Scales 

At distances r < 1 pc from the Galactic center, the gravitational potential is dom- 
inated by stars and by the supermassive black hole (SBH). The stars are observed 
to have a density^ 

p*(r) « 3.2 x lO 5 M pc- 3 (r/lpc)- 14 , r < 5pc. (4) 

Estimates of the SBH's mass range from 2 - 4 x lO 6 M p2EEEl the most recent es- 
timates, based on the orbits of single stars, suggest ~ 4 x 1O 6 M0, but for consistency 
with most pre-2005 papers on dark matter in the Galactic center, we here adopt 
Mbh = 3.0 x 10 6 M Q . The gravitational influence radius r^ of the SBH (the radius 
containing a mass in stars equal to twice the SBH's mass) is 77, » 2 pc. Within r^, 
the dark matter distribution is likely to have been strongly affected by whatever 
processes resulted in the formation of the SBH and the nuclear star cluster, and by 
any subsequent interactions between dark matter and stars. 

Fig. 1 shows a number of possible models for the dark matter density on sub- 
parsec scales. The standard halo model (SHM) discussed above, normalized to p = 
0.3 GeV cm -3 at the Solar circle, predicts p(rh) ~ lO 15 M pc _3 . Allowing for 
adiabatic contraction of the dark matter by the baryons (§2) increases this by a 
factor of ~ 10 2 , and the inner density slope 7 = — dlogp/dlogr in this contracted 
model is ~ 1.5 or steeper, implying a divergent annihilation flux. But even steeper 
density profiles are possible. If the SBH grew in the simplest possible way - via 
slow, symmetrical infall of gas - the density of matter around it would also grow, 
in the same way that contracting baryons steepen the overall Galactic dark matter 
scenario predicts a final density (of stars or dark matter) near the 

SBH of 

/ r \-7-j ! 
p(r)»p(r v ) — , 7sp = 2 + -— (5) 

\ r sp J * — 7c 

with 7 C the dark matter density slope in the pre-existing "cusp," i.e. the region 
r ^ i"h*^ Such a spike in the dark matter density would extend inward from 
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r sp « 0.2rftP^ and for j c > 1, 7 sp > 2.3 (Fig. 1). Now, the stellar density profile 
would respond in a similar way to growth of the SBH, and is known to be 

shallower than ~ r~ 23 (Eq. But this is probably a result of dynamical 
evolution that occurred after the SBH was in place (see below). 

At even smaller radii, a limit to the dark matter density is set by self- 
annihilations: p < p a = m/((Tv)t, with t w 10 10 yr the time since formation of the 
spike!^ A strict inner cut-off to the dark matter density is set by the black hole's 
event horizon, rs c h — IGM^jc 2 ss 3 x 10~ 7 pc for a non-rotating hole, although 
for reasonable values of m and (cry), the density is limited by self-annihilations well 
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Fig. 1. Possible models for the dark matter distribution in the Galaxy. The thin curve shows the 
standard halo model (SHM), and the thick curve is the same model after "adiabatic compression" 
by the Galactic baryons (stars and gas). Both curves are normalized to a dark matter density 
of 0.3 GeV cm -3 at the Solar circle. Curves labelled "spike" show the increase in density that 
would result from growth of the Galaxy's (SBH) at a fixed location. The annihilation plateau, 
P = pa = m/{av)t, was computed assuming m = 200 GeV, (av) = 10 — 28 cm 3 s _1 and t = 10 10 
yr. Dashed vertical lines indicate the SBH's Schwarzschild radius (rg c h m 2.9 X 10 -7 pc, assuming 
a mass of 3.0 x 10 6 M o and zero rotation), the SBH's gravitational influence radius (r^ ?s 1.7 pc), 
and the radius of the Solar circle (Rq ~ 8.0 kpc). Effects of the dynamical processes discussed in 
this article (scattering of dark matter off stars, loss of dark matter into the SBH, etc.) are excluded 
from this plot; these processes would generally act to decrease the dark matter density below what 
is shown here, particularly in the models with a "spike." 



We will discuss in the next section the prospects of observing DM annihilation 
radiation from the Galactic center. The dependence of the annihilation signal on 
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the DM profile is usually contained in the factor J defined as 



J An = KAfl- 1 / dip / p 2 dl (6) 

with Afi the angular acceptance of the detector, dl a distance increment along the 
line of sight, and ip the angle between the line of sight and the Galactic center; the 
normalizing factor K is normally set to K~ x = (8.5kpc)(0.3GeV/cm 3 ) 2 . Henceforth 
we denote by J3 and J5 the J values corresponding to Af2 = 10~ 3 sr and 10~ 5 sr 
respectively; the former is the approximate ang ular acceptance of EGRET^ while 
the latter corresponds approximately to atmospheric Cerenkov telescopes (ACTs) 
like VERITA^, CANGAROC03 and HES^El and to the proposed satellite ob- 
servatory GLAST.ES 

In dark matter models with an inner spike it is possible to solve analytically the 
integral of p 2 , and the resulting expression for J is 



J An 



3 An R% \ 2 7sp 

^ Jo P 2 j\ 1 
3 AQ R% 1 - 3/(2 7sp ) 

10 f Pa 





2 / \ 3 



, , (?) 

with r a the outer radius of the region where the density is limited by self anni- 
hilations (Fig. 1); these expessions assume that Aft ^S> (r sp /i?Q) 2 and the latter 
two expressions assume r a <§; r sp . A feeling for the range of plausible J values can 
be had by computing the p a and r a values corresponding to a set of ten, minimal 
supergravity benchmark models!^! Setting p Q — 0.3 GeV cm~ 3 and j sp — 2.4, and 
assuming that the dark matter cusp follows p oc r~ 1,5 outside of the spike and in- 
side the Solar circle (the "adiabatically compressed" version of the SHM, one finds 
1.2 x 10 4 < AilJ < 1.5 x 10 6 . In the absence of the spike, J values are many orders 
of magnitude lower. Such a wide range of J values makes it difficult to constrain m 
or (eru) from a measured annihilation flux or even to conclude that a signal would 
be detectable. 

Fortunately, many of the models in Fig. 1 can be ruled out. Once the dark matter 
distribution has been set up, it will evolve, due to interactions between dark matter 
particles, stars, and the SBH. In most (though not all) circumstances, this evolution 
has the effect of decreasing the dark matter density at r < . The most important 
evolutionary mechanisms are: 

• Scattering of dark matter particles off of stars. Stars in the Galactic nucleus 
have much larger kinetic energies than dark matter particles, and gravita- 
tional encounters between the two populations tend to drive them toward 
mutual equipartition, (l/2)m*v 2 « (l/2)mv 2 . Since the stars are the dom- 
inant population (at least at the current epoch), the dark matter heats up, 
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/ M bh \ 1/2 / r h \ 3/2 / to* , 



with to* the mean stellar mass. This heating tends to lower the dark mat- 
ter density and leads ultimately (in a time of several ?h eat ) to a density 
profile of the f orm p ~ r -3/2 < r < ^ flattens a pre-existing 

Thoat is also roughly the time for exchange of kinetic energy 
between stars, as a result of which the stellar density profile itself e volves 
toward a steady-state form, though with a steeper index p* ~ r - 7 / 4 El£5l 
This is probably the origin of the power-law stellar density cusp at the 
Galactic center (eq. 0}; the observed index, —1.4, is not quite as steep as 
the theoretical value but is probably consistent given the measurement un- 
certainties and given that the Galactic center contains stars with a range 
of masses and luminosities. 

Interestingly, even if the stellar and da rk matter cusps were once destroyed 
by the scouring effect of a binary SBH,E5E3E3I both might have been re- 
generated by this mechanism. In the case of the stars, as long as the density 
within the SBH's influence radius remains large enough that the star-star 
relaxation time is shorter th an ~ 10 10 yr, the r~ 7 / 4 cusp will re-form via 
the Bahcall-Wolf mechanism EH. This is likely to be the case if the mass 
ratio of the pre-existing binary SBH was extreme, e.g. ~ 10 : 1. Once 
the stellar cusp is back in place, heating of dark matter particles by stars 
will drive the dark matter towards its steady-state distribution, p ~ r~ 3 / 2 , 
although with a lower normalization than before cusp destruction. While 
many such evolutionary scenarios are possible (and should be worked out 
in more detail), the existence of a dense, collisional nucleus of stars like 
that observed in the Galactic center is strong circumstantial evidence of a 
steeply-rising dark matter density near the SBH. 

Capture of dark matter within the SBH. Any dark matter particles on orbits 
that intersect the SBH are lost in a single orbital period. Subsequently, 
scattering of dark matter particles off stars drives a continuous flux of dark 
matter into the SBHp2l Chan ges in orbital angular momentum dominate 
the flux; in a time Th ea t, most of the dark matter within will have been 
lost, although the net change in the dark matter density profile will be 
more modes t th an this suggests since more particles are continuously being 
scattered in!^H 

Capture of dark matter within stars. Another potential loss term for the 
dark matter is capture wit hin stars, due to scattering off nuclei followed by 
annihilation in stellar cores ! I I I However this process is not likely to be 
important unless the cross section for WIMP-on-proton scattering is very 
large. 



I 



I 



February 2, 2008 17:43 WSPC/INSTRUCTION FILE ms 



8 Bertone and Merritt 

These effects, as well as dark matter self-annihilations, can be modelled in a 
time-dependent way via the orbit-averaged Fokker-Planck equation!^ 



df 
dt 



1 dF E 
~4tt 2 p dE 



f{E)u coil {E) - f(E)u lc (E) 



(9) 



with f(E) the phase-space mass density of dark matter, E = —v 2 /2 + <f>(r) the 
energy per unit mass of a dark-matter particle, and <j)(r) the gravitational potential 
generated by the stars and SBH. The terms on the right-hand side of Eq. @ de- 
scribe the diffusion of dark matter particles in phase space space due to scattering 
off stars (Fe); loss of dark matter due to self-annihilations and/or capture within 
stars (Vcoll); and the loss-cone flux into the SBH {v\ c . Eq. © assumes an isotropic 
distribution of dark-matter velocities; this assumption is likely to break down very 
near the SBH, but the angular-momentum dependence of the loss cone flux is well 
understood and can be approximated via an energy-dependent loss term ^(E 1 )!^ 
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Fig. 2. (a) Evolution of the dark matter density at a radius of 10 — 5 r'h ?s 2x 10 pc from the center 
of the Milky Way due to the physical processes discussed in the text (Eq. |5J. Initial conditions 
were the standard halo model discussed in the text plus spike; the upper(lower) set of curves 
correspond to an initial density normalization at of 10(100) Mq pc — 3 . Mass and annihilation 
cross section of the dark matter particles were set to m = 200 GeV, (av) = 10 _27 cm 3 s _1 . In 
order of increasing thickness, the curves show the evolution of p in response to heating by stars; to 
self-annihilations; and to both processes acting together. Time is in units of T heat denned in the 
text; t = 10 corresponds roug hly to 10 10 yr. (b) Evolution of J averaged over an angular window 
of 10~ 5 sr. 



Fig. 2a shows the evolution of the dark matter density computed in this way, 
starting from a p ~ r~ 2 33 spike (p ~ r _1 cusp). Two values were taken for the 
initial density normalization at r = r^, p(fh) = (10, 100)M©pc~ 3 ; these values 
bracket the value p(rh) « 30M Q pc -3 obtained by extrapolating the SHM inwards 
from the Solar circle with p(R&) — 0.3GeVcm~ 3 (Fig. 1). The self-annihilation 
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term in Eq. @ was computed assuming m = 200 GeV, (av) = 10 _27 cm 3 s — The 
early evolution is dominated by self-annihilations but for t > 10 9 yr « Th ea t, heating 
of dark matter by stars dominates. The change in J(Ai7 = 10 -5 ) (Fig. 2b.) is 
dramatic, with final values in the range 10 4 < J < 10 5 . 



4. Direct and Indirect Detection 

In order to understand the nature of the dark matter, it is crucial to search for 
non- gravitational signatures. Physics beyond the standard model is actively being 
investigated using accelerators, and this is one of the main goals of the upcom- 
ing Large Hadron Collider (LHC), expected to begin operation around 2007 with 
proton-proton collisions at center-of-mass energies of 14 TeV. Numerous classes of 
models which provide interesting dark matter candidates will be tested at the LHC 
(Refe-ED-ESJ. 

An alternative approach is provided by so-called direct detection experiments. If 
the Galaxy is filled with WIMPs, then many of them should pass through the Earth, 
making it possible to look for the interaction of such particles with baryo nic : ma tter, 
e.g. by recording the recoil energy of nuclei as WIMPs scatter off themPEl The 
signal in this case depends on the density and velocity distribution of WIMPs in 
the solar neighborhood and on the WIMP-nucleon scattering cross section. WIMPs 
can couple either with the spin content of a nucleus, through axial-vector (spin de- 
pendent) interactions, or with the nuclear mass, through scalar (spin-independent) 
scattering. The second interaction typically dominates over spin-dependent scatter- 
ing in current experiments, which use heavy atoms as targets. More than 20 direct 
dark matter detection experiments are either now operating or are currentl y in de- 
velopmen t. P resently, the best direct detection limits come from the CDMS^and 
Edelweiss experiments, which probe nucleon-WIMP scattering cross sections of 
order 10 _7 pb, for a WIMP mass of order 100 GeV. These experiments have ruled 
out the WIMP discovery claimed by the DAMA collaboration, ' although it may 
still be possible to find exotic particle candidates and Galaxy halo mod els which 
are able to accommodate the data from all current experiments ! ' ^ ' "I i ( I 

Alternatively, dark matter can be searched for indirectly, via the study of the 
local flux of positrons and antiprotons. In fact, if the paradigm of WIMPs as massive 
particles in thermal equilibrium in the early universe is correct, dark matter particles 
are expected to annihilate in the Galactic halo producing possibly large fluxes of 
secondary particles. The High Energy Antimatter Telescope (HEAT) observed a 
flux of cosmic positrons well in excess of the predicted rate, peaking around ~ 
lOGeV and extending to higher energies!^ which could be the product of dark 
matter annihilations if sufficient clumpin g we re present in the halo fRefs.l7 T? | -l8ti[) . 
Upcoming experiments, such as AMS-02(S3 PAMELA,EHI and Bess Polar will 
refine the positron spectru m consid erably. An antiproton signal could also provide 
a signature of dark matter lgUgiJg£l 

Unlike other secondary particles, neutrinos produced by dark matter annihila- 
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tions can escape from dense media in which such annihilations may take place. For 
example, WIMPs that are captured in the Sun or Earth can annihilate at great 
rates. Although gamma rays cannot escape these objects, neutrinos often can, pro- 
viding an interesting signature to search for with high-energy neutrino telescopes 
fRefs. """3*1- 110111 . Although strongly model-dependent, the limits on the neutrino 
flux can be used to set constraints on dark matter particles. 

Dark matter could also be indirectly detected through its annihilation radiation 
in the Galactic halo or in extra-galactic sources. The prospects of detecting syn- 
chrotron radiation due to dark matter annihilations in the Galactic center have been 
studied in Refs. 11021 - 11041 while Ref. 11051 contains a discussion of the prospects for 
detecting the neutrino flux. Annihilation radiation could be enhanced by the pres- 
ence of substructures in the Galactic halo, either "clumps" (Refs. 11061 - I112|l or 
"caustics" (Refs. I"""j~3~l - 11150 . Recently, high-resolution numerical simulations have 
pointed to the existence of mini-hal os w ith masses as small as 10~6 solar masses 
and sizes as small as the solar system ESI (see also Refs. lll7llll8|l . If confirmed, this 
could have important consequences for indirect DM searches, and certainly deserves 
further investigation. 

The gamma-ray extra-galactic background produced by dark matter annihila- 
tions taking place in all structures and substructures in the Universe has been 
investigated in Refs. 11191 - 11211 One usually finds that the prospects to observe 
this signal are less promising than in the case of the gamma-rays from the Galactic 
center however (see e.g. Ref. I124|l . 

Much attention has been devoted to the study of gamma-rays from dark matter 
annihilations in the Galactic center fRefs. "T"""*"l - ll27[l . Given a detector with angular 
acceptance AO sr, the observed flux of photons produced by annihilations of dark 
matter particles of mass m and density p i d 1 

w-^S 7 - (io) 

where dN/dE is the spectrum of secondary photons per annihilation, (av) is the 
velocity-averaged self-annihilation cross section, dl is an element of length along the 
line of sight and Jao was defined in Eq. ("""J. The flux from solid angle Af2 is then 

<">(AQ, E) « 1.9 x 10- 12 f in - 26 M 3 r (^V^AOcm-Vi. (11) 
dE 10 *cm L \ m J 

Inward extrapolation of the density of the standard halo model (SHM) from the 
Solar circle into the Galactic center gives J 3 as 10 2 and J 5 w 10 3 . These J- values 
are large enough to produce observable signals for m any interesting choices of {av} 
and m. On the other hand, a recent detect ior P^ I by ACTs of a Galactic center 
gamma-ray source with energie s up to 1 TeV requires J values that are 10 2 to 10 4 
times bigger than in the SHMpEESIl assuming that the signal comes from dark 
matter annihilations. Achieving such large J values without abandoning the ACDM 
paradigm would require a significant enhancement in the dark matter density very 
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near the Galactic center, for instance, in the form of a density "spike" near the 
SBH. 

As discussed in §3, a steeply-rising dark matter density will due to self- 
annihilations and to dynamical interactions with the baryons. The results from 
a large set of time integrations of Eq. JHJ are summarized in Table 1 . Two extreme 
particle physics models were considered. In the first case, in order to maximize 
the ratio {av)/m, a cross section (cru)th = 3 x 10~ 26 cm 3 s -1 was assumed and 
m was set to 50 GeV. Higher values of (crv), though possible, would imply a low 
relic density, making the candidate a subdominant component of the dark matter 
in the universe. The lower limit on the mass strictly applies only to neu tralin os in 
theories with gaugino and sfermion mass unification at the GUT scale while 
the limit on the mass of KK particles is higher. The second extreme case assumed 
(<rv) = 0. Table 1 gives values of 7(Aft = 10~ 3 ) and J(Aft = 10~ 5 ) at t = 10 for 
each of these extreme particle physics models and for a variety of initial condtions. 
The final J-values depend appreciably on the particle physics model only when the 
initial dark matter density has a spike around the SBH; in other cases the central 
density is too low for annihilations to affect J. Particularly in the case of maximal 
(crv), the final J values are modest, log 10 ^3 ~ 5. 3 and log 10 J 5 < 7.0, compared 
with the much larger values at r = in the presence of spikes. 



Table 1. Results of time-integrations of Eq. J^J from a variety of initial dark 
matter models. 7 C and -y Bp are the power-law indices of the initial cusp and spike 
respectively. r c is the dark matter core radius in units of r^ f» 2 pc; p on r~ 1,/2 
for r sp < r < r c . Density at Rq is in units of GeV cm -3 . J 3 and J 5 are values 
of J averaged over windows of solid angle 10~ 3 sr and 10~ 5 sr respectively and 
normalized as described in the text. The final two columns give J in evolved models 
for (ov) = (no annihilations) and for ({crv). M) = (3 X 10 _26 cm 3 s _1 , 50GeV) 
(maximal annihilation rate), respectively. 















logio J 3 (Js) 




Model 


7c 


7s P 


r c 


p(R Q ) 


T — 


T = 10 


T — 10 


1 


1.0 






0.3 


2.56(3.51) 


2.56(3.50) 


2.56(3.50) 


2 


1.0 






0.5 


3.00(3.96) 


3.00(3.94) 


3.00(3.94) 


3 


1.0 




10 


0.3 


2.54(3.33) 


2.54(3.33) 


2.54(3.33) 


1 


1.0 




100 


0.3 


2.38(2.65) 


2.38(2.65) 


2.38(2.65) 


5 


1.0 


2.33 




0.3 


9.21(11.2) 


3.86(5.84) 


2.56(3.52) 


6 


1.0 


2.33 




0.5 


9.65(11.7) 


4.31(6.29) 


3.00(3.96) 


7 


1.0 


2.29 


10 


0.3 


6.98(8.98) 


2.61(3.88) 


2.54(3.33) 


8 


1.0 


2.29 


100 


0.3 


5.98(7.98) 


2.39(2.99) 


2.38(2.65) 


9 


1.5 






0.3 


5.36(7.30) 


4.81(6.58) 


4.78(6.53) 


10 


1.5 






0.5 


5.80(7.75) 


5.26(7.03) 


5.23(6.98) 


11 


1.5 




10 


0.3 


4.51(5.82) 


4.51(5.82) 


4.51(5.82) 


12 


1.5 




100 


0.3 


3.85(4.23) 


3.85(4.23) 


3.85(4.23) 


13 


1.5 


2.40 




0.3 


14.3(16.3) 


8.81(10.8) 


4.81(6.58) 


14 


1.5 


2.40 




0.5 


14.8(16.8) 


9.25(11.3) 


5.25(7.02) 


15 


1.5 


2.29 


10 


0.3 


9.67(11.7) 


4.77(6.51) 


4.51(5.82) 


1(5 


1.5 


2.29 


100 


0.3 


7.67(9.67) 


3.87(4.64) 


3.86(4.23) 
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Table 2. Parameters in the fitting function for the boost. 







AU = 10" 3 








AU = 10~ 5 






Model 


-Brain 


Bmax 


a 


6 


-Bmin 


Bmax 


a 


b 


5 


-0.02 


1.31 


0.66 


0.73 


-0.05 


2.35 


0.55 


1.50 


6 


-0.01 


1.31 


0.66 


0.51 


-0.06 


2.34 


0.56 


1.28 


7 


-0.02 


0.05 


0.75 


0.92 


-0.18 


0.38 


0.72 


1.31 


8 


-0.18 


-0.17 


0.76 


1.37 


-0.86 


-0.51 


0.73 


1.64 


13 


2.14 


6.28 


0.42 


-0.05 


2.95 


7.35 


0.40 


0.10 


11 


2.16 


6.29 


0.43 


-0.28 


2.97 


7.36 


0.41 


0.13 


15 


1.96 


2.21 


0.74 


-0.27 


2.32 


3.02 


0.71 


0.07 


16 


1.30 


1.31 


0.75 


0.53 


0.73 


1.14 


0.73 


0.84 



In these evolutionary models, the value of J at a fixed time t — rThcat is deter- 
mined completely by the initial conditions and by the quantity (crv) jm that specifies 
the annihilation rate. This outcome can be expressed in terms of the boost factor b 
defined as J/Jjv, with Jjv the value in the SHM hav ing t he same density normal- 
ization at r = Rq as in the evolving model. One find J ^ ' that the boost factors at 
r = 10 (roughly 10 10 yr) can be well approximated by the function 

B(X) = B max - {l/2)(B max - B mm ){\ + tanh[ a (A - b)}} (12) 

with X = log 10 ((CTv)/10- 30 cm 3 s- 1 )/(TO/100GeV) and B = log 10 b. Values of the 
fitting parameters are given in Table 2 for the "spike" models of Table 1; as noted 
above, in the absence of a spike, the final J values are essentially unaffected by 
annihilations hence inde pendent of (crv)/m. Recent 

analyseJEBESD f the HESS 
Galactic center dat JP^suggest that the observed gamma-ray spectrum is consistent 
with particle masses of order 10 TeV and cross sections of order 3 x 10 _26 cm 3 s _1 , 
if boost factors are as high as 10 3 < b < 10 4 . Table 2 shows that such boost factors 
are indeed achievable at X = log 10 (3 x 10" 26 /10 _30 )/(10 4 /100) w 2.5 if the initial 
dark matter distribution is similar to that of Models 13 or 14, i.e. a "spike" inside 
of a p ~ r" 15 cusp. 

Annihilation radiation might also be detected from the cent ers o f galaxies other 
than the Milky Way fRefs. IT3*3*1 - I143[) . Even globular clusters^EI have been sug- 
gested as possible targets. As in the case of the Milky Way, the major uncertainty 
in predictions of the annihilation flux from external galaxies is the unknown distri- 
bution of dark matter on small scales; absent power-law cusps or spikes, indirect 
detection of dark matter from ext ernal galaxies is probably impossible with the 
current generation of detectors As argued above, a plausible guide to the dark 
matter distribution on small scales is the stellar density profile. The Local Group 
galaxies M31 and M32 both exhibit steep stellar cusps, ~ r ~2* , 1 -5 < 7+ < 2.0, 
inward of ~ 1 pc, similar to what is observed in the Milky Wayj a l these galaxies 
also show clear kinematical evidence for SBHs. By contrast, the Local Group galax- 
ies M33 and NGC205 exhibit cores in the s tellar density and any SBHs in these 
galaxies are too small to be detected J2SE3 Stellar cusps like the one in the Milky 
Way would be unresolved at the distance of the Virgo cluster, but it is a reasonable 
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guess that all comparably-luminous galaxies contain similar distributions of lumi- 
nous and dark matter on sub-parsec scales. The situation is likely to be different for 
giant galaxies like M87 in Virgo, which are known to have much shallower stellar 
density profiles within their SBH's sphere o f infl uence, probably a consequence of 
an earlier epoch of scouring by binary SBHsP^ The matter distribution near the 
centers of giant galaxies is probably essentially unchanged since the last round of 
galaxy mergers, making them unfelicitous sites for indirect detection of dark matter. 
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